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ABSTRACT 

Observations of redshifted 21-cm radiation from neutral hydrogen during the epoch of reion- 
ization (EoR) are considered to constitute the most promising tool to probe that epoch. One of 
the major goals of the first generation of low frequency radio telescopes is to measure the 3D 
21-cm power spectrum. However, the 21-cm signal could evolve substantially along the line 
of sight (LOS) direction of an observed 3D volume, since the received signal from different 
planes transverse to the LOS originated from different look-back times and could therefore 
be statistically different. Using numerical simulations we investigate this so-called light cone 
effect on the spherically averaged 3D 21-cm power spectrum. For this version of the power 
spectrum, we find that the effect mostly 'averages out' and observe a smaller change in the 
power spectrum compared to the amount of evolution in the mean 21-cm signal and its rms 
variations along the LOS direction. Nevertheless, changes up to ~ 50% at large scales are 
possible. In general the power is enhanced/suppressed at large/small scales when the effect is 
included. The cross-over mode below/above which the power is enhanced/suppressed moves 
toward larger scales as reionization proceeds. When considering the 3D power spectrum we 
find it to be anisotropic at the late stages of reionization and on large scales. The effect is 
dominated by the evolution of the ionized fraction of hydrogen during reionization and in- 
cluding peculiar velocities hardly changes these conclusions. We present simple analytical 
models which explain qualitatively all the features we see in the simulations. 

Key words: cosmology: theory, cosmology: diffuse radiation, cosmology: reionization, meth- 
ods: numerical, methods: statistical 



1 INTRODUCTION 

The epoch of reionization (EoR), when the first luminous sources 
reionized the neutral hydrogen in the intergalactic medium (IGM), 
is currently the frontier of observational astronomy. Observations 
of the Cosmic Microwave Background Radiation (CMBR) (Ko- 
matsu et al. 2011; Larson et al. 2011) and high redshift quasar ab- 
sorption spectra (Becker et al. 2001; Fan et al. 2006; Willott et al. 
2009) jointly suggest that reionization took place over an extended 
period spanning the redshift range 6 ^ 2; ^ 15 (see e.g. Mitra et al. 
2011). Observations of high redshift Ly a-emitting galaxies (Mal- 
hotra & Rhoads 2004; Ouchi et al. 2010; Kashikawa et al. 2011) 
and gamma ray bursts (Totani et al. 2006) are also consistent with 
this picture. 

Observations of redshifted 21-cm radiation are considered to 
constitute the most promising tool to probe the EoR (for a review 
see Furlanetto et al. 2006). For the past few years substantial efforts 
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have been undertaken both on the theoretical and experimental side 
(reviewed in Morales & Wyithe 2010). The first generation of low 
frequency radio telescopes ( GMRT\ LOFAR^, MWA^, PAPER^) 
is either operational or will be operational very soon. Preliminary 
results from these facilities include foreground measurements at 
EoR frequencies (Ali et al. 2008; Bernardi et al. 2009; Pen et al. 
2009; Paciga et al. 2010) as well as some constraints on reioniza- 
tion (Bowman & Rogers 2010; Paciga et al. 2010). 

Motivated by the detection possibility of the EoR 21-cm signal 
and the subsequent science results, a wide range of efforts are ongo- 
ing on the theoretical side with the goal to understand the physics of 
reionization and its expected 21-cm signal. Furlanetto et al. (2004) 
developed analytical models to calculate the ionized bubble size 
distribution and use this for calculating the 21-cm power spectrum. 



^ http://gmrt.ncra.tifr.res.in/ 

http://www.lofar.org/ 
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Such models are very useful in predicting the signal quickly for a 
wide range of scales and investigating the large parameter space. 
However, they cannot incorporate details of reionization and be- 
come less accurate when the bubbles start overlapping. Numerical 
simulations are probably the best way to predict the expected 21- 
cm signal. Although challenging, there has been considerable pro- 
gresses in simulating the large scale 21 -cm signal during the entire 
EoR (Iliev et al. 2006; Mellema et al. 2006a; McQuinn et al. 2007; 
Shin et al. 2008; Back et al. 2009). More approximate but much 
faster semi-numerical simulations of the structure and evolution of 
reionization and the 21 -cm signal have also been developed (Zahn 
et al. 2007; Mesinger & Furlanetto 2007; Santos et al. 2008; Geil 
& Wyithe 2008; Thomas et al. 2009; Choudhury et al 2009). These 
methods are capable of generating volumes with sizes as large as 
~ IGpc^ (Alvarez et al. 2009; Santos et al. 2010). Many aspects 
such as source properties, feed back effects, distribution and prop- 
erties of sinks have also been investigated in detail (see Trac & 
Gnedin 2009, for a review on reionization simulations). 

One of the major goals of all first generation EoR telescopes 
is to measure the spherically averaged 3D 21 -cm power spectrum. 
Measurements of the 21 -cm power spectrum will provide a wealth 
of information about the timing and duration of reionization, large 
scale distribution of H I and its evolution, source properties and 
clustering (AH et al. 2005; Sethi 2005; Datta et al. 2007; Lidz et al. 
2008; Barkana 2009). To obtain the spherically averaged 3D power 
spectrum one needs to average over the 3D volume produced by the 
observations. Of this 3D volume one axis (the LOS axis) is along 
the frequency direction. Since light from the lower frequency side 
of the 3D volume takes a longer time to reach us than light from 
the high frequency side, the observer will see reionization in an ear- 
lier phase at the lower frequency side than at the higher frequency 
side. The statistics of 21 -cm fluctuations could therefore be chang- 
ing over the observed volume. As we will see in section 2 and 3, 
in some reionization scenarios the change could be substantial es- 
pecially near the end of reionization. Almost all previous studies 
calculate the 3D 21 -cm power spectrum without taking this effect 
into account. In this paper we investigate the effect of LOS evo- 
lution or the so called 'light cone' effect on the measured 21 -cm 
power spectrum, using numerical simulations to quantify it. Under- 
standing the light cone effect is important because it will be present 
in the data and needs to be taken into account when interpreting the 
observed 21 -cm power spectrum. Our aim is to understand under 
which conditions and at what scales this effect needs to be consid- 
ered. 

The light cone effect is well known from studies of galaxy 
clustering (see e.g. Matsubara et al. 1997). In the context of 21-cm 
studies of reionization it was first considered by Barkana & Loeb 
(2006). These authors studied analytically the anisotropic structure 
of the two point correlation function caused by the effect. This 
appears to be only work that considered the effect of a changing 
source population. However, more work has been done on the light 
cone effect for a single bright source, such as a QSO. For this case 
the effect will make the H II region appear to be teardrop-shaped 
(Wyithe et al. 2005; Yu 2005; Majumdar et al. 2010). The effects 
on the power spectrum and correlation function for this case were 
investigated by Sethi & Haiman (2008). In addition, for very lu- 
minous sources the effect of relativistically expanding H II regions 
(Shapiro et al. 2006) would have to be added to the one purely due 
to evolution of the signal along the LOS. 

Bright QSOs are quite rare, so the more common form of the 
light cone effect will be due to the evolving source population and 
the growing H II regions around groups of sources. Our aim is to 



study this version of the effect on the spherically averaged 3D and 
the ID LOS power spectra using reahstic numerical simulations of 
reionization. 

The paper is organized as follows. Section 2 briefly describes 
our simulations and the procedure used to generate light cone 
cubes. We present our results in Section 3. Section 4 describes 
two simple toy models which explain qualitatively the main fea- 
tures we see in the simulation results. Section 5 investigates how 
the inclusion of peculiar velocities affect our results. We summa- 
rize our results and conclusions in Section 6. The cosmological 
parameters we use throughout the paper are Qrn = 0.27, Q^k = 
0, Qb = 0.044, h = 0.7, n = 0.96 and erg = 0.8, consistent with 
the WMAP seven-year results (Komatsu et al. 201 1). 



2 SIMULATION 

2.1 The redshifted 21-cm signal 

The 21-cm radiation is emitted when neutral hydrogen atoms go 
through spin-flip transitions. The radiation can be decoupled from 
the cosmic microwave background (CMB) photons either through 
collisions with hydrogen atoms and free electrons (Purcell & Field 
1956; Field 1959; Zygelman 2005) or through hya photon pump- 
ing (Wouthuysen 1952; Field 1959; Chen & Miralda-Escude 2004; 
Hirata 2006; Chuzhoy & Shapiro 2006). This makes 21-cm radia- 
tion detectable either in emission or absorption against the CMB. 
The differential brightness temperature with respect to the CMB is 
commonly written using the spin temperature Ts as 

5T.«27.4.H.mK (^^y^Jl^^^il + Sn) (1) 

where xm and are the mass averaged neutral fraction and the 
density fluctuations of hydrogen. Note that the 21-cm signal re- 
mains undetectable when the spin temperature Ts is coupled to the 
CMB temperature Tcmb- During the EoR, Ts is expected to be 
coupled to the gas kinetic temperature through hya photon cou- 
phng. In addition the gas kinetic temperature is expected to be 
much higher than the CMB temperature due to heating by shocks. 
X-rays and Lya photons. This would make the redshifted 21-cm 
signal visible in emission. We assume here that Ts ~ Tgas ^ 
Tcmb which makes the 21-cm signal independent of the actual 
value of Ts. This is a reasonable assumption during the later stages 
of the EoR. During the initial stages of reionization, when there are 
few sources of radiation, this assumption might not hold (Back et 
al. 2010; Thomas & Zaroubi 2011). The next subsection describes 
how we simulate the fluctuations in the H I density. 



2.2 N-Body and radiative transfer runs 

Details about our simulation methodology (N-body simulation and 
the subsequent radiative transfer) have been presented in previous 
papers (lUev et al. 2006; MeUema et al. 2006a; lUev et al. 2007, 
2011). Iliev et al. (2011) described the simulations we use here in 
more detail. Here we only present a brief overview of the major 
features of these simulations. 

We start by simulating the evolution of the dark matter distri- 
bution using the CubeP'^M N-body code in a comoving volume of 
(163cMpc)^ using 3072^ particles and 6144^ cells. This implies 
particle masses of 5.5 x lO^M© and a minimum resolvable halo 
mass of ~ lO^M© which approximately matches the minimum 
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Figure 1. The evolution of the mean mass averaged neutral fraction xm and the rms of 21 -cm fluctuations with redshift for the two different reionization 
simulations LI and L3. 




k Mpc"^ 



k Mpc 



Figure 2. The evolution of dimensionless spherically averaged 3D 21 -cm power spectrum A'^j^(k) with redshift in the LI (left panel) and L3 (right panel) 
simulations. The light cone effect is not taken into account. 



mass of haloes able to cool by atomic cooling. The N-body simu- 
lations give the DM density field, locations and masses of haloes. 
We then assume that the baryons trace the DM density field and 
assign an ionizing photon luminosity to each halo assuming it to be 
proportional to the halo mass, 



Mh^b 



(2) 



where is the number of ionizing photons emitted per time. Mh 
and nip are the halo mass and proton mass respectively. The effi- 
ciency parameter can be written as. 



9i = 



lOMyr 



(3) 



where fj is the number of ionizing photons emitted into the IGM 
per baryon per star forming episode (which is taken to be the same 
as simulation time-step, about 11.5 x 10^ years). This makes the 
factor fj the product of the escape fraction of ionizing radiation 
/esc, the star formation efficiency /*, and the number of ioniz- 
ing photons produced per baryon for a given initial mass function 
(IMF), Nj. The latter number is around 5,000 for a Salpeter IMF 
and ^10 times higher for top heavy IMFs, and the two fractions 
are of the order 10%. This gives fj values in the range 1 to 100. 
We divide the halos into low mass atomically cooling halos of mass 



10^-10^ M0 and high mass atomically cooling halos of masses 
higher than that. To take into account feedback on the low mass 
halos, we turn off their ionizing luminosity when they are located 
in an ionized region. 

We then calculate the transfer of ionizing photons with the C^- 
Ray code (Mellema et al. 2006b) on a 256^ grid. Ionizing photons 
are traced from every source cell to every grid cell within a given 
time step. This gives us the distribution of the H I fraction in the 
volume at different redshifts. 

Here we consider three cases for the ionizing photon luminos- 
ity. In the first simulation, labeled LI in Iliev et al. (201 1), g^ — 8.7 
or = 10 for sources of mass > 10^ M© and g^ — 130 or 
= 150 for sources of mass between 10^ — 10^ M©. In the sec- 
ond simulation, labeled L2, these factors are g^ = 1.8 or = 2 



and g^ 



.7 or — 10. In the third case, labelled L3, sources 



ization at almost the same time as LI, the active sources have been 
assigned a higher luminosity with = 21.7 or = 25. This 
particular setup was chosen to make it analogous to the older sim- 
ulations from Iliev et al. (2008), but updated for the WMAP5 cos- 
mology. 

The simulations LI and L3 are the ones which have the 
strongest evolution, where L3 due to the lack of low mass sources 
has the fastest evolution. Since the light cone effect is caused by 
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Figure 3. Two-dimensional slices of the 21 -cm signal from a coeval cube at z = 9.31 (left) and light cone cube (right) with central redshift Zc = 9.31 
(simulation LI). For the latter the x-axis corresponds to the LOS direction. Comparison of these two slices shows the effect of evolution on the sizes of H II 
regions at both the front and back sides of the Hght cone cube. Note that for this visuaHzation the mean signal has not been subtracted from the light cone cube. 



evolution, below we will focus on these two cases and briefly men- 
tion the more slowly evolving case L2 at the end of Section 3. 

2.3 Light cone cube 

The simulations provide us with so-called 'coeval cubes', 3D vol- 
umes of density and H I fraction at the same cosmological red- 
shift. The extent of these cubes corresponds to a redshift range 
of Az ^ 0.6 — 0.9, depending on redshift. An observer can not 
observe these coeval cubes, but we use them to create observable 
'light cone' cubes. The procedure which was previously introduced 
in Mellema et al. (2006a), is as follows: 

• From the simulation we obtain a set of N coeval 21 -cm cubes 
at redshifts zi, Z2 zn{zi < zi < 2; at) each of integer size 

and physical comoving size . For the case at hand M — 256 
and L = 163 cMpc. 

• Starting at z\ we create a redshift series 2;lc of length m — 
K X M (K ^ N) which will constitute the redshift (LOS) axis 
of the light cone 'cube' (which will therefore not be cubical). Each 
consecutive redshift in the series is the same comoving distance 
apart, namely L/M. 

• We then construct the light cone cube by stepping through this 
redshift series and constructing the 21 -cm slices of size for 
each redshift. 

• To create the pth 21 -cm slice of the light cone cube, we first 
calculate the integer division p/M and its remainder q. We pick 
up the ^th slice from the two coeval cubes at zi and z^+i, where 
zi ^ zi,c{p) ^ zi-^i, and use linear interpolation in redshift to 
create an 21-cm slice at 2;lc(p)- 

We should point out that the light cone cubes constructed this 
way differ from the observational ones in that the field of view has 
a constant comoving size and not a constant angular size. This is 
a natural consequence of the way they are constructed from the 
simulation results and makes it easier to construct the 3D power 



spectra from them. For the real interferometric observations the an- 
gular field of view would be slowly changing as a function of fre- 
quency and the physical comoving size depends on the redshift via 
the angular- size distance relationship. For determining the 3D 21- 
cm power spectrum in /c-space it will always be possible to extract 
a volume with a constant comoving field of view from the observa- 
tional data. 

From the LI simulation we extracted 35 coeval cubes at red- 
shifts spanning from z = 11.20 to 8.4. In this model reionization 
starts earlier with the mass weighted ionization fraction reaching 
1% and 50% around redshifts z = 17.22 and 9.46 respectively 
(see Figure 1). In the L3 model, where the smaller mass halos do 
not contribute, the reionization process starts later (because mas- 
sive sources form later) and the 1% and 50% points are reached 
around redshifts z = 12 and 9 respectively. By construction, the 
two simulations complete reionization at the same redshift of 8.4, 
so in the L3 simulation reionization proceeds faster. Because L3 
has fewer sources, the characteristic bubble size for a given neutral 
fraction xhi is bigger. Both models are consistent with the recent 
CMB measurements of the electron scattering optical depth. Fig- 
ure 1 shows the evolution of the mass averaged neutral fraction xhi 
(left panel) and the rms of 21-cm fluctuations (right panel) for the 
two models. Note that for L3 the rms is higher than for LI because 
of the larger ionized bubbles which amplify the rms signal. 

Since the comoving distance between z — 11.20 to 8.4 is 
larger than 163 cMpc our full light cone cube is constructed by us- 
ing the periodicity of our cosmological volume. However, this does 
mean that we pass through the same structures several times and our 
power spectra would be unphysical below scales of ^ 0.08 Mpc~^ . 
We therefore limit our power spectrum analysis to subvolumes of 
LOS size 163 cMpc, which roughly corresponds to a frequency 
depth of - 10 MHz. 

Figure 2 shows the dimensionless spherically averaged 3D 21- 
cm power spectra Alj^{k) = A;^P3D(A;)/27r^ (Peacock 1999) for 
coeval cubes at different redshifts for both simulations. In the be- 
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Figure 4. The light cone effect on the 21-cm power spectrum A^^^ (k) in simulation LI. Left panels: A^^^ (k) for the light cone cube (red), coeval cube (blue) 
centered around the redshift mentioned in the plot and coeval cubes (two dashed lines) for redshifts corresponding to back and front sides. Right panels: the 
relative difference (A^^^^^ — Ag^^^^j/Ag^^^^ where 'cc' and 'Ic' stand for coeval and light cone cube respectively. 



ginning of reionization the power spectrum is dominated by the 
density fluctuations and quite similar to the underlying dark matter 
power spectrum. As reionization proceeds the growing ionized bub- 
bles add power at larger scales. When reionization reaches ^50% 
the power reaches a maximum at larger scales. As the neutral frac- 
tion goes down further the overall amplitude of the power spec- 
trum also goes down. Note once again that the L3 simulation has 
more power than the LI model because of the larger ionized bub- 
bles. From Figure 2 we see the power spectrum evolve both in am- 
plitude and in slope (see Lidz et al. 2008, for a detailed discus- 
sion on the power spectrum evolution). The details of the evolu- 
tion depend on the reionization scenario (for example, the evolu- 
tion is much faster in the L3 model). If we consider LI, we can see 
that the power spectrum Al-Q{k) at /c = 0.1 Mpc~^ changes from 
- 0.25 mK^ to - 10 mK^ in the redshift range z = 8.46 to 8.89. 
Such rapid evolution of the power spectrum (by a factor of ^ 40 at 
/c = 0.1 Mpc~^ ) within Az = 0.43 provided the motivation for 
studying the light cone effect. 



3 EFFECT OF EVOLUTION 

3.1 Spherically averaged power spectrum 

In this section we present and discuss our results on the light cone 
effect. Figure 3 shows two 21-cm images constructed from simula- 
tion LI, the left one from a coeval cube at z = 9.31, and the right 
one from a light cone cube, where the horizontal axis corresponds 
to the LOS direction and the central redshift is z = 9.31. We see 



Table 1. Details about the simulated cubes from simulation LI, used for our 
analysis. 



Zc 


a^Hlc 


Redshift 




rms 






extent 


variation 


variation 










(mK) 


8.76 


0.18 


8.46 - 9.07 


7 X 10-3 -0.33 


1.0 - 9.7 


9.31 


0.44 


8.99 - 9.65 


0.3-0.55 


9.5 - 10 


9.94 


0.63 


9.59 - 10.3 


0.54-0.68 


10 - 9.4 


11.2 


0.76 


10.8 - 11.63 


0.73-0.8 


9.3 - 9.2 



that ionized regions (black patches) at the higher z side (righthand 
side) are smaller in the light cone image than in the coeval image. 
Conversely, the ionized regions at the lower z end (lefthand side) 
are larger in the light cone cube. 

Figures 4 (for LI) and 5 (for L3) show the effect of evolution 
on the spherically averaged 3D 21-cm power spectrum. Note that 
we do not include the modes k(/ccc = 0, /c^ = 0, kz) when we take 
spherical average over all modes between k and k + dk. Since we 
are performing this analysis in the context of radio-interferometric 
measurements that do not measure the modes 'k{kx = 0, /c^ = 
0, it is appropriate not to include those modes. It is in fact 
quite important for the analysis we present. Excluding the modes 
k(A:x = 0, /c^ = 0, kz) changes the 3D power spectrum from the 
light cone cube considerably. We discuss this more extensively in 
Appendix A. 
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Details about the central redshift {zc), mass averaged neutral 
fraction at the central redshift (xhic), redshift range, neutral frac- 
tion range and the range of the rms variations in the 21 -cm signal 
for the different light cone cubes from simulations LI and L3 are 
given in Tables 1 and 2 respectively. The left panels of Figure. 4 
show the spherically averaged 3D power spectra /S.%j)(k) for light 
cone cubes at four different central redshifts. For comparison it also 
shows the power spectra for coeval cubes at three redshifts (high, 
low redshift end and central redshift of the light cone cube). The 
righthand panels plot the relative difference in the power spectra 



between the coeval cube at the central redshift and the light cone 
cube. General features we see in all panels are that the effect is 
stronger on large scales and increases as we go up in scale. In ad- 
dition, we find that the power is enhanced (suppressed) with re- 
spect to the coeval cube at large (small) scales for neutral fractions 
xm ^0.5. During the last phase of reionization (bottom panels) 
the power in the light cone cube is suppressed at all scales. At red- 
shift z = 11.20 we see that the power spectrum is enhanced by 
^ 15% at modes k < 0.5Mpc~^ and suppressed by ^ 10% 
at small scales. At redshift z = 9.94 we do not see much effect 
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Figure 7. Anisotropy of the 3D power spectrum. The plots show the ratio A^^^^^/ A|j^^^ as a function of jj? for different A:-modes for two central redshifts 
Zc = 8.76 (left panel) and 10.02 (right panel) for L3 simulation. 




except on the largest scale where the power is larger by 15%. 
This is rather surprising since the evolution of the 21 -cm signal is 
stronger in this redshift interval than in the z = 11.20 band. For 
the cube centered around redshift z — 9.94 the neutral fraction and 
the rms change from 0.547 to 0.68 and 10 to 9.4 mK and the power 
spectrum is amplified by a factor of ^ 7 at /c = O.lMpc"^ (see 
Table 1 and the last two left panels of Figure 4), much more than 
what we see in the cube around z = 11.20. So we would expect a 
larger effect than what we see in Figure 4 at z = 9.94. This trend 
continues and we see almost no effect for redshift z — 9.31 where 



the neutral fraction and the rms change even more (0.325 0.547 
and 9.4 10 mK) and the power spectrum is amplified by a factor 
of ^ 5 at /c = 0.1 Mpc~^ At redshift z = 8.76, instead of an 
enhancement we see suppression on all scales with differences up 
to 30%. In the L3 simulation this suppression is up to 50% at red- 
shift z — 8.76. All other features are quite similar in the L3 model 
(Figure 5) even though the reionization process proceeds faster and 
the ionized regions are larger in this model. 

Another way to describe the trend we see is that we find a 
cross-over mode /ccross-over below which power is enhanced and 
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Figure 9. Cartoon illustration of Toy model 1 . The left panel represents the coeval cube containing ionized bubbles of the same size. The right panel represents 
the light cone cube with bubbles from back/front side appearing smaller/larger in comparison to the coeval cube. 



Table 2. Same as Table 1 for simulation L3. 





a^Hic 


Redshift 




rms 






extent 


varies from 


varies from 










(mK) 


8.76 


0.27 


8.46 - 9.07 


1 X 10-2 - 0.5 


2.5 - 13.5 


9.31 


0.63 


8.99-9.65 


0.45 - 0.76 


13.2 - 13.3 


10.02 


0.85 


9.59 - 10.3 


0.75-0.9 


13.3 - 12.2 


10.57 


0.93 


10.2 - 10.96 


0.88-0.96 


12.4-11.7 



above which it is suppressed. The cross-over scales /ccross-over 
shifts towards lower k as the reionization proceeds. At the end of 
reionization the cross-over mode is lower than the lowest mode we 
measure from the simulation box. 

3.2 Light cone effect as a function of LOS width 

Above we present results using the entire cubes i.e, for a LOS width 
corresponding to the size of our simulation volume. However, it is 
interesting to explore how the effect changes as one reduces the 
LOS width. Obviously in the limit of small widths, the light cone 
effect will disappear, so considering a range a widths allows us to 
study how it varies with width. 

Here we consider sub-boxes of different LOS widths A 2; and 
calculate the quantity (Ag^cc — ^3dic)/^3d1c for different k 
modes. Figure 6 shows (AId^c — ^3dic)/^3d1c a function of 
LOS width (and Az) for different k modes at two central redshifts 
Zc — 8.76 and 10.02 for the L3 simulation. As expected, we see 
that the quantity (AId^c — ^3dic)/^3d1c decreases for smaller 
LOS width (and /\z). As discussed later in the subsection 4.3, we 
expect the quantity to increase quadratically with the LOS width. 
However, due to the smaller number of modes available for the 
smaller Az, the results are too noisy to test this expectation, al- 
though they are roughly consistent with it. 

We do not show results for the other two central redshifts of 
L3 and LI simulation as the light cone effect is relatively smaller 
for these, but find similar results there. These results suggest that 
measurements of the light cone effect for different LOS widths can. 



in principle, be used to correct for the effect or at least find the sign 
of the effect. 



3.3 Anisotropics in the power spectrum 

The light cone effect introduces an anisotropy in the full 3D 21 
cm power spectrum. For a fixed /c-mode, the power spectrum 
will depend on A;||, the LOS component of k. Peculiar velocities 
and the Alcock-Paczynski effect are the other major sources of 
anisotropics in the 21 cm power spectrum. In order to understand 
the anisotropic power spectrum and to separate the physics from 
astrophysics (Barkana & Loeb 2005) each effect should be studied 
in detail. Though the first generation of low frequency radio tele- 
scopes (i.g, LOFAR, GMRT, MWA ) are unlikely to able to mea- 
sure the anisotropics in the 21 cm power spectrum, this will be the 
ultimate goal of such measurements. 

Fig. 7 plots the ratio Ag^ic/AsDcc a function of ji^ for 
different /c-modes for two central redshifts Zc — 8.76 (left panel) 
and 10.02 (right panel) for L3 simulations. Here /i = k\\/k. In 
the left panel {zc — 8.76) we see that the ratio Ag^ic/AlDcc 
mcreases from - 0.7 (at /x^ = 0.1) to 1.1 (at /x^ = 0.9) for 
k — 0.16 Mpc~^. For k — 0.3Mpc~^, the ratio increases from 
^ 0.8 to 1 for the same ji^ range. For higher /c-modes the degree 
of anisotropy decreases and the power spectrum is becoming more 
isotropic. We do not see any significant anisotropics for the cen- 
tral redshift Zc = 10.02 (right panel) where the neutral fraction 
xhic = 0.86. The other redshifts of the L3 simulation also do not 
show significant anisotropics and the results for the LI simulation 
are similar to L3. 

We do not try to quantify the anisotropics further as we see 
the curves are not very smooth due to the small number of modes 
at large k. Our results are sample variance limited and should be 
considered as qualitative rather than quantitative. Larger simulation 
volumes are needed to quantify the anisotropics more precisely. 
Barkana & Loeb (2006) (fig. 2) reported significant anisotropics 
at large scales (r = 100 Mpc) for neutral fraction xhi = 0.25 
in their Pop III reionization model. This is qualitatively consistent 
with our results. 
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3.4 ID LOS power spectrum 

Recent results by Marker et al. (2010) show that the dimension- 
less ID LOS power spectrum A^^ lqs (^) ^ can be extracted more 
accurately as there is no large scale bias (which may arise due to 
foreground subtraction) and smaller error bars in the recovered ID 
LOS power spectrum (shown in figure 11 of Harker et al. 2010). 
The ID LOS power spectrum can also extend to smaller scales be- 
cause of the higher resolution in the frequency direction. Motivated 
by this, we also study the effect of evolution on the ID LOS power 
spectrum. Figure 8 shows the effect in the ID LOS power spectrum 
for the LI case. Interestingly, the ID LOS power spectrum is hardly 
affected by light cone effects, except near the end of reionization. 
We do not show results for the L3 simulation as these are very sim- 
ilar to LI. This result adds one more advantage to the strategy of 
measuring ID LOS power spectrum. However, the ID LOS power 
spectrum is very flat on large scales due to the aliasing of high k- 
modes. Therefore, although it can be more accurately measured, it 
may be difficult to see the impact of ionized bubbles and extract the 
reionization physics. 

We would like to mention here that besides the simulations LI, 
L3 we also analysed the L2 simulation (for more details see Iliev et 
al. 2011) where the reionization is much more gradual and overlap 
happens at a redshift z = 6.5. The evolution is thus relatively slow 
and we see that the light cone effect is smaller in amplitude but 
otherwise has similar features as what we presented above for the 
cases LI and L3. 



3.5 Comparison to previous work 

Barkana & Loeb (2006) analytically studied the light cone effect 
using the two point correlation function ^ (r, /i), rather than consid- 
ering the power spectrum. Since they considered different reioniza- 
tions models, redshift range and the cosmological parameters, as 
well as another diagnostic, we can here only provide a qualitative 
comparison. 

Barkana & Loeb (2006) limited their investigation to the late 
stages of reionization (xni < 0.5) and find that the effect is sig- 
nificant from the time when the reionization is ~ 70% complete 
to its end (see their fig. 4) and that large scales are affected more 
than small scales. We find similar results as we see the largest dif- 
ferences in the power spectra at large scales for the later stages of 
reionization. We also find substantial differences (10 — 30%) in 
the first half of reionization, but this phase was not considered in 
Barkana & Loeb (2006). 

For a fixed correlation length r Barkana & Loeb (2006) 
showed that the correlation function ^(r, ji) for )Lx = 1 is identi- 
cal to the value f or fi = around xm = 0.5, suppressed close to 
the end of reionization, and enhanced in the intermediate period. 
As explained above, /x = 1 corresponds to the LOS direction and 
thus measures the light cone effect. Therefore, for a fixed length 
scale r Barkana & Loeb (2006) found the light cone effect to have 
a negligible impact before and around xm = 0.5, to suppress the 
correlation function at the end of reionization but to enhance it in 
the intermediate period. This is exactly what we find. In the L3 sim- 
ulations, we find that the power spectra are suppressed during the 
late stages of reionization but enhanced before that. 



^ ^m,LOs(^) = ^^iD,LOs(fc)/7r (Peacock 1999) 



4 TOY MODELS 

To understand the results from the previous section we consider 
here two analytical toy models of reionization. 



4.1 Toy model 1 

In this toy model we consider a very simple scenario. We consider 
N number of spherical, non-overlapping and randomly placed ion- 
ized bubbles in a uniform H I medium in the coeval cube. The 
spherically averaged 3D power spectrum for such a scenario can 
be written as 

N 

PMk) = Y.V?W\kRi) (4) 

i=l 

where Vi — ^TvRi and W{kR) is the spherical top hat window 
function defined as 

W(kR) = — ^ [sm( kR) - kR cos( kR)] (5) 

k-^R-^ 

Now P3r>{k) = J2iLi for k < l/i?max, where i?max is the 
radius of the biggest bubble in the cube since W{x) ^ 1 for x < 1. 

For the coeval cube we assume that all bubbles are of the same 
size Vo, therefore the power spectrum can simply be written as, 

P3ncc{k) = NV^. (6) 

Because of the evolution effect in the light cone cube, bubbles at 
the back side will appear smaller and bubbles at the front side will 
appear bigger and in addition their shapes could be somewhat elon- 
gated along the LOS (see Figure 1 of Majumdar et al. 2010). To 
make our calculations simpler we assume the bubbles in the light 
cone cube are spherical but have different sizes Vo + AVi. As we 
saw in the simulation results, the global ionization fraction for light 
cone cubes is almost the same as in the coeval cube at the central 
redshift, so we assume X^^^ AVj = 0. The spherically averaged 
3D power spectrum for the light cone cube at larger scales is then 
given as, 

N 

PsDicW = ^(K + Ay,)' 

i=i 

N N 

= PsDcc (k) + 2Vo JZ^^'^Yl 

i=l i=l 

N 

= P3Dcc{k) + Y,^y^ (7) 
i=l 

The above equation explains two major features we see in the simu- 
lation. First it explains why the light cone effect is relatively small. 
We see that the effect cancels out in the linear order. Only the 2nd 
order term J2iLi^^^ survives the averaging and affects the light 
cone power spectrum. So in this sense the light cone effect is a 
'2nd order effect' in the spherically averaged power spectrum. Sec- 
ond, because J2iLi always positive the power is always 
enhanced at larger scales which is exactly what we see in the sim- 
ulation. When the bubble sizes are not identical in the coeval cube, 
we can subdivide the entire range of bubble sizes into small size 
bins. The above analysis would then be applicable to each individ- 
ual size bin and thus to the entire range. 

Our second toy model considers a slightly more realistic but 
still quite simple reionization scenario. 
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kMpc 

Figure 10. The relative change in the 3D 21-cm power spectrum /S!^j^{k) 
due to the combined effect of linear growth of structure and adiabatic ex- 
pansion of the Universe. The ordinate is the quantity 100 x (A^^^^ — 
^3Dwo)/^3Dwo' where 'w/wo' means the effect is/is not included. We 
use very high redshift simulation cubes (before the reionization starts) to 
calculate this. 
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Figure 11. Bubble size distribution calculated from the simulated cubes 
around redshift z = 9.09 for L3 model. The dashed lines show the bubble 
distributions for three coeval cubes at redshifts z = 9.38, 9.09 and 8.76 
(from left to right) corresponding to the back, middle and the front side of 
the light cone cube centered around redshift 9.09. The solid line (red) show 
the distribution for this light cone cube. 



4.2 Toy model 2 

In this toy model, we consider a reionization model in which spher- 
ical ionized bubbles of different sizes are placed randomly. If there 
is no overlap between ionized bubbles then the ionized fraction 
would be, 



/ 



(8) 



where n{R) is the number density of bubbles of size R. But in 
practice randomly placed bubbles will overlap with each other and 
expand further to conserve the emitted photon numbers. We neglect 
the fact of further expansion of bubbles and so the actual ionized 
fraction which would be less than the above can be calculated using 
(Furlanetto et al. 2004) 



Xi^l - exp(-(5). 



(9) 



The spherically averaged 3D power spectrum can be calculated 
from the two point correlation function using the relation 



P3D(fc) =47r r drr^arf'''^'''^^ 

Jr=0 



kr 



(10) 



Here we would like to mention that the evolution makes the correla- 
tion function ^ anisotropic i.e, a function of both r and /x (Barkana 
& Loeb 2006; Sethi & Haiman 2008). Since our aim is to study the 
light cone effect on the spherically averaged power spectrum and 
qualitatively understand the main features we see in the simulation 
results, we assume ^ to be isotropic. We expect this approximation 
not to affect our conclusions. The correlation function ^(r) can be 
decomposed into (Zaldarriaga et al. 2004) 



(11) 



where ^xx, ^ss and ^xS are the correlation functions of ionization 
field, density field and the cross-correlation between two fields re- 
spectively. 

In the density field correlation function ^ss, two quantities 
change with redshift: 1) density fluctuations grow with time, 2) the 
mean density decreases because of the expansion of the Universe. 
Thus the evolving ^ss in principle would contribute to the light 
cone effect. In the linear regime the two quantities together scale 



as ^ (1 + z)~^'^. For a distance of 163 cMpc they jointly change 
^ 3% along the LOS. We use very high redshift simulation cubes 
(before the reionization starts) and find < 0.1% enhancement in 
the 3D power spectrum on almost all scales (Figure 10). This re- 
sult agrees with McQuinn et. al (2006) who predicted a constant 
enhancement in the power spectrum (see their Appendix A). The 
contribution of the evolving ^ss to the total light cone effect on the 
power spectrum is therefore negligible and hence we ignore this 
term in the rest of our analysis. 

The evolution of ^ is thus mainly dominated by the ^xx on 
scales larger or comparable to the size of ionized bubbles. For the 
rest of our analysis we only consider the term ^xx- The function 
which is defined as ^xx{r) = {X1X2) — should be zero both for 
= and Xi = 1. It also should satisfy the boundary conditions 
(for details see Zaldarriaga et al. 2004) 



Xi 





for r 
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The correlation function can be calculated for a given bubble 
distribution as (Furlanetto et al. 2004) 



{xiX2){r) = 

1 — exp 



+ exp 
X I 1- 



/ 



exp 



m'^v.iR) 



. j dR^[V(R)-Vom 



(12) 



where Vo{R, r) is the volume of the overlap region between two 
ionized regions centered a distance r apart. The function can be 
written as 



Vo{R,'> 



A-kR^/2, - nr[R^ - r^/12] for r < 2R 
for r > 2R 



In the next subsection we present some bubble size distribu- 
tions measured from simulation. We will then model the bubble 
distribution and use that for the subsequent analysis. 
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Coeval Cube 
Light Cone Cube 



Table 3. Parameters for the bubble size distributions in Toy model 2. 
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Figure 12. Bubble size distributions used to calculate the power spectrum 
in the Toy model 2. SoHd (red) and dashed (green) lines represent coeval 
and light cone cubes respectively. 
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Figure 13. The power spectra for the different bubble size distributions 
from Figure 12. The y-axis is scaled arbitrarily. 



4.2.1 Bubble size distribution and its evolution 

We calculate the bubble size distribution from the simulation using 
the spherical average method (see Friedrich et al. 2011, for more 
details on different bubble size estimates). Fig. 1 1 shows the bubble 
size distribution RdP/dR^ for coeval cubes at three redshifts cor- 
responding to the back, middle and front side of a light cone cube 
centered around redshift 9.09. It also shows the bubble size dis- 
tribution in the light cone cube centered around 9.09. The coeval 
distributions at the three redshifts differ considerably. For example 
the radii at which the bubble distribution peaks are 10, 20, and 90 
cMpc. Interestingly the bubble size distribution for the light cone 
cube is very similar to the coeval box at the central redshift. In the 
light cone cube the bubble size distribution would be the average 
of those of the coeval cubes in the redshift range Zc =b /S.z where 
2 is the extent of the cube along redshift axis. Although the bub- 
bles in the light cone cube are smaller/larger in the back/front side 
compared to the coeval cube, the average bubble distribution in the 
whole light cone cube is very similar to the coeval cube of the cen- 
tral redshift. Because of this 'averaging effect' the light cone effect 

^ This is quantity is essentially same as V{R){dn/dR) and follows the 
condition Xi = f (dP/dR) dR 
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is small even though there is a substantial evolution in the bubble 
distribution across the box. 

We investigate further and to make the following calculations 
simpler we parameterize V^R)-^^ — Aexp( — (i? — Rc)^ /aji). 
This is motivated by Fig. 1 1 (see also Figure 2 in Furlanetto et al. 
2004). We normalize the function using Eq. 8. 

Now consider the case around the central redshift Zc. Since 
the light cone cube covers the redshift range Zc ± it will have 
slightly more bubbles both at the large and small bubble size ends 
than the coeval cube at redshift Zc. This is exactly what we see in 
the simulation (Figs. 3 and 11). 

As we mentioned in Section 4.1 the average neutral fraction 
xhi for a coeval cube at redshift Zc and for a light cone cube cen- 
tered around redshift Zc is almost the same. We consider reioniza- 
tion for two values of Q = 0.5 and 0.7 (see Eq. 8). Parameters for 
the bubble distribution are summarized in Table 3. Figure 12 shows 
the bubble distribution for (J = 0.5 (solid) and 0.7 (dashed) for the 
coeval cube and the light cone cube. As we discussed above, there 
will be more large and small size bubbles in the light cone cube 
compared to the coeval cube, we approximate this by increasing 
an for the light cone cube. The bubble size at which the quantity 
V-^^ peaks has been kept same for both for a fixed Q. We also 
see in the Fig. 12 that for (J = 0.5 the number density of bubbles 
of size Rb > IS cMpc is higher in the light cone cube than the 
coeval cube. The 'cross over radius' i.e, the bubble radius beyond 
which the number density becomes higher than in the coeval cube 
is IScMpc. For Q = 0.7 the cross over radius 27cMpc) is 
higher than for (J = 0.5. This is because for higher Q the char- 
acteristic bubble size increases. Obviously, the exact distribution 
could be different but the general features such as the increase of 
characteristic bubble size and cross over radius for larger Q values, 
and larger bubbles in the light cone cube than the coeval cube, are 
likely to be true in all reionization scenarios where stars/QSOs are 
dominant sources. Since our aim is to qualitatively understand the 
effect of evolution, we can use these simplified distributions. 

Figure 13 plots the power spectrum for the bubble distribution 
models we describe above. We see that there is a scale /ccross-over 
below (k < /ccross-over) wMch thc light cone cube has more 
power than the coeval cube and above (k > /ccross-over) which 
it is the other way around. This is because the number of bub- 
bles in the light cone cube below the cross-over radius is less 
than in the coeval cube. We call this scale the 'cross-over mode'. 
We see this feature in Fig. 13 where the cross-over modes are 
fcross-over = 0.097 Mpc"^ and 0.063 Mpc"^ for Q = 0.5 and 
0.7 respectively. The cross-over mode is thus seen to shift towards 
larger scales (smaller k) as Q increases. This is because the cross- 
over radius is larger for (J = 0.7. Based on these results we find the 
empirical relation between the cross-over radius and the cross-over 
mode to be 

1.7 

-Rcross — over — ~. (1-^) 



"Across — over 
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This simple toy model explains the two main features seen in 
the simulation results 

(i) The power spectrum in the light cone cube is en- 
hanced/suppressed on large/small scales compared to the one from 
the coeval cube at the central redshift. 

(ii) The cross-over mode shifts towards large scales as reioniza- 
tion proceeds. 

As we pointed out in Sect. 3, we do not see a large effect when the 
mean ionization fraction is around 50%, even though the evolution 
is across the light cone cube is substantial at that stage. We can now 
understand this to be because the cross-over mode shifts towards 
larger scales as reionization proceeds and around 50% ionization 
the cross over scale is already almost the same as lowest mode that 
we can measure from our simulation volume, even though it has a 
size of 163 cMpc. We therefore predict that for a larger simulation 
volume, enhanced power on scales k < 0.08 Mpc~^ will be found. 

4.3 Taylor expansion of the power spectrum evolution 

In addition to the more heuristic models given above, the follow- 
ing approach also helps in understanding some of the trends we 
see in Figure 4 and 5. We find that the coeval cube power spec- 
trum A|e)cc(^5 z) for a given mode k changes very smoothly with 
redshift z. So we expand Ago^cl^, ^) ^ Taylor series around a 
central redshift Zc as 

ALcc(/c,^) = ALcc(/c,^c)+a(AL)+6(AL)'+c(AL)'+..,(14) 
where AL is the comoving distance from the central redshift 

. » - - p— . = ^ - 

h (^^?^), ' ^ = (^W^)^ • We ignore the higher 
order terms. Next we calculate the light cone power spectrum 
A3Dic(^5 Zc) by taking average of Al^ccC^, z) in the range zbL/2 
around redshift Zc using 

ALic(/c,z) = y ALcc(/c,z)^/L, (15) 

^ J-L/2 

where L is the comoving LOS width. The above equation can be 
simplified to 

ALic(/c,z) = ALcc(/c,Zc) + 6^. (16) 

We see that the linear term (aAL) and all terms with odd pow- 
ers cancel out and only the quadratic term (&(AL)^) and the other 
terms with even powers survive the averaging process. This sup- 
ports our argument that the light cone effect is a '2nd order effect' 
and that linear trends in the evolution of the power spectrum av- 
erage out. The fractional change in the power spectrum due to the 
light cone effect is given by — . Positive/negative values of the 

^^^3Dlc 

parameter b denote that light cone power is suppressed/enhanced 
compared to the coeval value. 

To test this quadratic approximation we use simulation L3 and 
fit the polynomial Eq. 14 for a given mode k around three different 
central redshifts, taking L = 163 cMpc. Using the values of b we 
calculate the percentage change in the power spectrum and also 
measure the actual percentages from the simulation results. Values 
for three different k modes are given in Table 4. From these it can 
be seen that the quadratic expansion correctly predicts the sign of 
the parameter b and reproduces the trends seen in the simulations. 
During the early phases the match with the simulations is quite 



Table 4. Comparison between the quadratic expansion and simulation re- 
sults. Listed are the relative sizes of the light cone effect as predicted by the 
quadratic expansion (Q) and measured in the simulation (S). 
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good, but at later stages it under predicts the changes. Most likely 
the discrepancies are due to the neglect of the higher order terms. 



5 EFFECT OF PECULIAR VELOCITY ON THE LIGHT 
CONE EFFECT 

The peculiar velocity of the IGM gas influences the 21 -cm power 
spectrum (Bharadwaj et al. 2001; Bharadwaj & Ah 2004). During 
the dark ages when the H I density is expected to trace the DM 
density, the spherical averaged power spectrum is enhanced by a 
factor of 1.87 at linear scales. As reionization proceeds, the relative 
contribution of the peculiar velocity to the 21 -cm power spectrum 
changes considerably with redshift. For inside-out reionization sce- 
nario the peculiar velocity could increase the 21 -cm power spec- 
trum by a factor of 5 (see Mao et al. 201 1, Fig. 3) during a short 
period in the beginning of reionization (xi < 0.2, see Mao et al. 
201 1, Fig. 3). When reionization is at its ~50% phase, peculiar ve- 
locity effects slightly decreases the 21 -cm power spectrum on the 
large scales relevant for the first generation of EOR experiments. In 
other words, peculiar velocity effects change the evolution of 21- 
cm power spectrum and hence could affect the light cone effect. We 
briefly investigate this here. 

The method for taking the peculiar velocity into account when 
constructing the light cone cube was outlined in Mellema et al. 
(2006a) and described in detail in Mao et al. (2011); in the ter- 
minology of the latter we use the MM-RMM(1 xRT) scheme. The 
left panel of Figure 14 shows the 21 -cm power spectrum with pe- 
culiar velocity for coeval cubes at three different redshifts (center 
and two ends) as well as for the light cone cube. The right panel 
shows the relative difference (A^dcc — ^suic) / ^suic The figure 
looks mostly very similar to Figure 4 where we did not include any 
peculiar velocity effects, the exception being the earliest stages, at 
redshift z = 11.20. Here the case with peculiar velocity shows 
a negligible light cone effect whereas the case for no peculiar ve- 
locity shows a ~ 10% difference in the power spectrum at large 
scales. The reason for this is that in the beginning of reionization 
bubble growth is relatively slow and evolution is dominated by the 
peculiar velocity. As reionization proceeds the evolution is mainly 
dominated by the growth of ionized bubbles and hence the peculiar 
velocity has almost no impact on the hght cone effect. Note that this 
does not mean that the inclusion of peculiar velocity does not affect 
the power spectrum. In fact peculiar velocity causes comparable or 
even larger changes than the light cone effect. A more thorough 
exploration of the effects of peculiar velocity will be presented in 
future work. 
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Figure 14. The effect of peculiar velocity (pv) on the light cone (Ic) effect in the LI case at different redshifts. Left panels: Dot dashed, dashed and solid 
lines show the power spectrum for cubes without pv but with Ic, with pv but without Ic and with both pv and Ic, respectively. The right panels plot ( A|j^^ 
^3Dic) /^3Dic where we incorporated the peculiar velocity in both A| 



^Idcc and Aj^^^. 



6 CONCLUSIONS AND DISCUSSION 

We investigate the effect of evolution on the 3D and ID LOS 21 -cm 
power spectra during the entire period of reionization. We use three 
different EoR simulations in a volume of 163 cMpc on each side, 
one in which reionization is more gradual, ending at z ~ 6.5 and 
two in which it is more rapid, ending at z ^ 8.5. In one of these 
rapid simulations, reionization is driven by more massive sources, 
leading to relatively larger ionized bubbles. Below we summarize 
our results: 

• For the cases we studied, the spherically averaged power spec- 
trum changes up to ^ 50% in the k range 0.08 to 9 Mpc~^ using a 
redshift interval corresponding to the full extent of our simulation 
volume, (163 cMpc). As expected, for smaller redshift bins the ef- 
fect is found to be smaller. Large scales are affected more and the 
effects at smaller scales are minor. 

• Substantial evolution of the mean mass averaged neutral frac- 
tion xhi, rms variations in the 21 -cm signal, and bubble size distri- 
bution along the LOS axis are averaged out in the spherically av- 
eraged power spectrum. This averaging effect makes the light cone 
effect relatively small compared to the evolutionary changes along 
the LOS axis. 

• We can detect anisotropics in the the full 3D power spectra 
on large scales in the later stages of reionization, but are unable 
to quantify the /i-dependence of this effect with the simulations 
available to us. 

• The bubble size distribution in the light cone cube centered 
around redshifts Zc is remarkably similar to the bubble size distri- 
bution in the coeval cube at Zc, even if there is substantial evolution 
in the ionized fractions along the LOS. This is the reason why we 
see a relatively small effect on the 21 -cm power spectrum com- 
pared to the amount of change in the xhi and the rms of the 21 -cm 
signal. 

• The large scale power is enhanced and the small scale power is 
suppressed most of the time except at the final phase of reionization 



where the power spectrum is suppressed at all scales we can mea- 
sure in our simulations. In other words, there is a 'cross-over mode' 
A^cross-over bclow and above which the power is enhanced and 
suppressed respectively. The cross-over mode /ccross-over shifts to- 
wards lower /c-mode (large scale) as reionization proceeds. 

• Surprisingly we see very little effect when reionization is 
~ 50% complete and there is a rapid evolution in the xhi and 
the rms. We argue that at this stage of reionization the cross-over 
mode /ccross-over is already comparable to the lowest k mode we 
can measure from the simulation and enhancement of power should 
be present at larger scales than that. 

• Despite the fact that the reionization histories differ consider- 
ably between the three simulations, we see quite similar results. 

• Growth of structures with redshift and the expanding back- 
ground enhance the power spectrum by ^ 0.1% for our 163 cMpc 
cube. Its evolution is therefore dominated by the ionization field 
during the reionization. 

• An analytical toy model (Toy model 1) can explain the large 
scale power enhancement due to light cone effect as well as its 
smallness. 

• A second analytical toy model (Toy model 2) for a light 
cone cube with more large bubbles beyond some cross-over radius 
i?cross-over and Icss bubblcs below that, can explain all the features 
we see in the simulation results. 

• The presence of more large bubbles and fewer small bubbles 
of size < i?cross-over is rcsponsiblc for the enhanced/suppressed 
power on scales below/above /ccross-over- The fact that the cross- 
over scale shifts towards lower k as reionization proceeds is be- 
cause the cross-over bubble size i?cross-over increases as reioniza- 
tion proceeds. 

• Interestingly we find that the light cone effect is less prominent 
in the ID LOS power spectra. 

We should note that instruments such as LOFAR and MWA 
are expected to measure down to k ^ 0.01 Mpc~^, scales larger 
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than we were able to analyze here (/cmin = 0.08 Mpc~^). From our 
results we expect enhanced power on those scales in the light cone 
cube. Especially when reionization is around ~ 50% we expect 
more enhanced power on these larger scales. Reionization simula- 
tions of even larger cosmological volumes would be useful to better 
understand the effects at those scales. On the other hand, the afore- 
mentioned telescopes will not reach beyond k ^ 1 Mpc~^ making 
the small scale light cone effects observationally less relevant. 

The removal of the large foreground signals of the EoR 21cm 
signal is expected to affect the large scale LOS modes k\\ signifi- 
cantly. Although details about which scales will be affected depend 
on the subtraction technique used, it is obvious that if L is the co- 
moving length over which the foreground subtraction is performed, 
modes with k\\ < 27r/L cannot be extracted (McQuinn et. al 2006). 
The equivalent bandwidth for the simulation boxes we consider is 
~ 10 MHz and it is likely that foreground subtraction techniques 
will use considerably larger bandwidths (see e.g. Chapman et al. 
2012). The same authors also show that foreground residuals do not 
affect the extraction of the 3D spherically averaged power spectrum 
over bandwidths of 8 MHz. However, the effects of foregrounds re- 
main clearly an issue which requires careful consideration when 
considering LOS effects in the 21cm signal. 

In our simulations the spherically averaged power spectra are 
based on equal numbers of modes in the LOS and transverse di- 
rections. However, most of the ongoing and upcoming surveys will 
not sample the full range in the spatial and frequency directions for 
many /c -modes. This is due to the fact that they have better resolu- 
tion in the frequency (LOS) direction than in the spatial directions. 
For example the LOFAR core has a maximum baseline which cor- 
responds a maximum transverse mode /c_Lmax ^ lMpc~^. The 
intrinsic frequency resolution of the array is better than 1 kHz, but 
likely the observed data will be stored with ~ 10 KHz frequency 
resolution, equivalent to LOS mode k\\ ^ 35Mpc~^. When 

^ II max ^ 

using this resolution to calculate the spherically averaged power 
spectra, the LOS modes k\\ < k will contribute more compared to 
the transverse modes for k > A:^niax- Since the light cone effect 
makes the power spectra anisotropic i.e, different power for differ- 
ent combintions of {k±,k\\) for a given k, the lack of small scale 
transverse modes k± ~ A;, in principle, would affect the power 
spectra measurements at those /c -modes. However, as shown in Fig. 
7, small scales are hardly anisotropic due to the light cone effect so 
we do not expect those modes to be affected much due to the in- 
complete sampling of small scale modes. In addition, small scales 
k > 0.6 Mpc~^ are expected to be dominated by system noise and 
are unlikely to be measured. 

Based on our results, we conclude that the light cone effect 
is important especially at scales where the first generation of low 
frequency instruments are sensitive. It can bias cosmological and 
astrophysical interpretations unless this effect is understood and in- 
corporated properly. 
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Figure Al. Illustration of the effects of including k(kx = 0, ky = 0, kz) 
modes on the 3D power spectrum (simulation LI). The abscissa shows the 
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and dot-dot-dashed lines are for redshifts 8.76, 9.31, 9.94 and 1 1.20 respec- 
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APPENDIX A: EFFECT OF INCLUSION OF 

k(kx = 0, ky = 0, kz) MODES ON THE POWER 
SPECTRUM 

Radio interferometric experiments cannot measure the modes at 



kx 



where k 



2'k(u,v) 



u,v are two compo- 



nents of the baseline vector U and r is the comoving distance. 
In order to predict the expected 21 -cm power spectra for some 
reionization model or interpret the observed 21 -cm power spec- 
tra the modes k(kx = 0, ky = 0, kz) should be excluded when 
the power spectrum is calculated from the simulated data. Fig- 
ure Al plots 100 X [PsB-inik) - P3D-ex(A;)]/P3D-ex(A:) with 
k for different redshifts for light cone cubes. Here P3D-in(A^) and 
P3D-ex(A:) are the 3D power spectra including and excluding the 
k(kx = 0, ky = 0, kz) modes respectively. We find that power is 
enhanced by 10 - 200% for A: < 0.1 Mpc"\ The reason is that for 
the light cone cube there is a gradual change in the mean bright- 
ness 6Tb with redshift. Large scale LOS modes with kx — ky — ^ 
gain power because of this and hence affect the large k modes in 
the spherically averaged power spectrum. We find that the coeval 
cubes are hardly affected because the mean 5Th is similar for all 
slices. We also note that for the simulations studied in this paper, 
exclusion of the k(kx = 0, ky = 0, kz) modes is practically the 
same as the subtraction of the mean brightness temperature from 
each single frequency 21 -cm map. 



